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При расчетах активной биологической защиты представляется полезным получить уравнения движе- 
нии релятивистских протонов в веществе при наличии магнитного поля, поскольку комбинированные 
способы защиты могут оказаться оптимальными в тех случаях, когда вес защитных экранов должен быть 
ограничен сверху. 

Как известно, уравнение движения заряженной частицы в магнитном поле при отсутствии вещест- 
ва в векторной форме имеет вид: 
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где т - масса частицы; е - заряд; с - скорость света в вакууме; ѵ - вектор скорости; В - вектор магнит- 
ной индукции. 


Поскольку нас будут интересовать частицы, скорость которых соизмерима со скоростью света, то 
уравнение (1) следует записать в форме соотношения 
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где р = Ш - импульс частицы, масса которой зависит от энергии, т. е. скорости частицы. 

При прохождении частицы через вещество будут иметь место различные взаимодействия частицы со 
средой, в результате чего энергия частицы будет уменьшаться вдоль траектории движения. В интересую- 
щем нас диапазоне энергий протонов главными являются потери энергии на ионизацию атомов защит- 
ной среды. Потери энергии на ионизации будут эквивалентны действию некоторой силы, действующей 
на частицу (протон) в направлении, противоположном направлению вектора скорости. Величина этой си- 
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лы будет равна удельной энергии ионизации на элемент длины траектории /, т.е. сила равна — . 
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При учете силы сопротивления вещества за счет ионизации уравнение (2) может быть представлено в 
векторной форме: 
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где прямые скобки означают абсолютное значение соответствующей величины, а ѵ - это модуль вектора 
скорости частицы. 

В релятивистском случае левая часть уравнения (3) должна быть представлена следующим образом: 
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где /п 0 - масса покоя частицы. 

В трехмерной прямоугольной системе координат имеет место соотношение 

ѵ 2 = х 1 + у 2 + (5) 

где для сокращения записи первые производные по времени обозначены точками над соответствующими 
координатами. 
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Проекция векторного уравнения (4) на оси координат с учетом (5) дает возможность записать уравне- 
ние движения в следующем виде: 
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Здесь В„ В у и В, - составляющие вектора магнитной индукции по осям координат. 
Для сокращения записи введем обозначение 
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и правые части уравнений системы (6) представим в виде 
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Раскрывая производные в левых частях системы (6) и учитывая (8), получим 
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Здесь р = - = 
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- релятивистский фактор, а две точки означает производную второго по- 
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так как все координаты являются функциями времени. Подставляя ( 10) в систему уравнений (9) и выпол- 
няя несложные преобразования, получим 
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Здесь введены для сокращения записи обозначения: 
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При моделировании обыкновенных дифференциальных уравнений на аналоговых вычислительных 
машинах (АВМ) принято разрешать эти уравнения относительно старших производных (метод пониже- 
ния порядка производной). 

Рассматриваем систему (11), как систему трех уравнений с тремя неизвестными производными 
Вычисляем главный определитель системы (11) 
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Применяя правило Саррюса или какой-либо другой способ вычисления определителей и учитывая (5) 
в несложных, но громоздких преобразованиях, получим 


А = с 6 (1-/3 2 ) 2 . (14) 


Теперь выражение, например, 

длях 

будет иметь вид 
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Громоздкие алгебраические преобразования при вычислении (15) позволяют в конечном итоге полу- 
чить сравнительно простое выражение для х 
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Аналогично вычисляются выражения для у и г. В результате получаем систему уравнений движения в 
виде: 
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т,'х = (-1 В г у - ± В..г | ]/Г=Т 2 - 9 (г-) хѴЦ- ?)*-> 
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При /3 = 0 получаем уравнения движения заряженной частицы в веществе при наличии магнитного по- 
ля для нерелятивистского случая 
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Полученные уравнения движения имеют форму, удобную для моделирования на АВМ. В цилиндри- 
ческой системе координат уравнения движения в релятивистском случае будут иметь вид: 
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При моделировании на АВМ полезно также учитывать связи между кинетической энергией частицы и 
ее скоростью в виде 


ѵ_ = Кб(б — 2б 0 ) ■ 
с б + б 0 




Здесь б 0 = ШцС 2 - энергия покоя частицы. 
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